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ABSTRACT 

We present the first survey of electric field data using the ARTEMIS spacecraft in the solar wind to study 
inertial range turbulence. It was found that the average perpendicular spectral index of the electric field depends 
on the frame of measurement. In the spacecraft frame it is -5/3, which matches the magnetic field due to the 
large solar wind speed in Lorentz transformation. In the mean solar wind frame, the electric field is primarily 
due to the perpendicular velocity fluctuations and has a spectral index slightly shallower than -3/2, which is 
close to the scaling of the velocity. These results are an independent confirmation of the difference in scaling 
between the velocity and magnetic field, which is not currently well understood. The spectral index of the 
compressive fluctuations was also measured and found to be close to -5/3, suggesting that they are not only 
passive to the velocity but may also interact nonlinearly with the magnetic field. 

Subject headings: magnetic fields - magnetohydrodynamics (MHD) - plasmas - solar wind - turbulence 



1. INTRODUCTION 

The solar wind is a plasma that is observed to be 
turbulent with fluct u ations at a broad range of scales 
(iTu& Marschl 11995b iGoldstein et aD 11995b iHorburv et al l 
120051: iBruno & Carbonel 120051: iPetrosvan et al.ll2010t) . It is 
usuall y modeled as a cascade of energy from large scales 
(e.g. JWicks et alj2010t) : where the energy is injected, to small 
scales (e.g. JChen et alj|2010al). where kinetic pr ocesses dis- 
sipate the energy (e.g., Schekochihin et al. 2009). The iner- 
tial range fluctuations are thought to be prim arily Alfvenic in 
natur e, with Alfven-wave-like polarizations (Belcher & Davis l 
1971 1 ) and phase speeds close to the Alfven speed (Ba le et al] 
2005). 

There are various theories of Alfvenic turbulence, based 
on interacting packe t s of A lfven waves. The theory of 
lGoldreich~& Sridhar (119951) . based on critical balance, 
predicts that the Alfvenic fluctuations have a perpendic- 

-5/3 

ular one-dimensional energy spectrum E(k^) ~ kj_ , 
where k± is the wave vector perpendi c ular t o the mag- 
netic field. The theory of iBoldyrevl J2006), which in 
addition assumes scale-dependent alignment, predicts 

-3/2 

that their spectrum is E(k±J ~ k ± . Similar predictions 
also exist for the multitude of imbalanced theories (e.g., 
Lithwick et al. 2007; Beresnvak & Lazarian 2008; Chandran 
20081: iPerez & Boldvrevl 120091: iPodesta & Bhattacharieel 
20ldlPodestall2011h . 

In the solar wind at 1 AU, it has been shown that the spec- 
tral index of the magnetic field is close to -5/3 on average but 
that the spectral index of the velocity is c l oser to -3 /2 (e.g. , 
Mangenev et al. 2001; Podesta et al. 2007; Salem et al. 200' 



Tessei n et al.l 12009; Pod esta & Borovsk v 2010; Wic ks et al 
2011). This difference between the two fields is not consis 



tent with any of the current theories of Alfvenic turbulence 
and is one of the currently unsolved problems of solar wind 
turbulence. 

Past measurements of the electric field spectrum in the 
frame of the spacecraft found it to closely m atch the magnetic 
field dBale et alj|2005b ISahraoui et al.ll2009l) . These measure- 
ments used single i ntervals of data but it has been shown (e.g., 
iTessein et al.l2009l) that the velocity and magnetic field have a 
large spread of spectral indices and many intervals are needed 



to determine the average behavior. 

In this Letter, we present a survey of electric field measure- 
ments in the solar wind using many intervals of data. We 
explain why the electric field in the spacecraft frame follows 
the magnetic field and make new measurements of the electric 
field in the mean solar wind frame. In Section [2] we describe 
the data set, in Section[3]we discuss our results and in Section 
2]we present our conclusions. 

2. DATA SET 



We used data from the ARTEMIS mission ( Angelopoulos 
1201 Oh. which is an extension of the THEMIS mission 
(Angelopoulos 2008). During late 2010, the two ARTEMIS 
spacecraft (PI and P2) moved from equatorial Earth orbits to 
Lunar Lagrange orbits (~ 60 Re from the Earth). Periods of 
solar wind data were selected in which each spacecraft was 
upstream of the Moon, out of Earth's ion foreshock and the 
required instruments were operational. The selected days are: 
days 245-257, 308-310, 316-318, 337-343 of 2010 and days 
1_3, 40-42 of 201 1 forP7; days 217-230, 275-284, 304-307, 
361-364 of 2010 and days 25-28 of 2011 for P2. The same 
day in both spacecraft was avoided so that the intervals are 
independent. All of the data from these days were split into 
6 hr sections resulting in 272 intervals, 98% of which were in 
slow solar wind (<500 km s" 1 ). 

Spin resolution (^3 s) ele ctric field data, E SI . , from the 
electric field instrument (EFI; Bonnell et al. 2008) was used, 
along with spin resolution magne tic field data, B, from 
the fluxgate magnetometer (FGM; lAuster et al.l 120081) and 
varying resolution ion velocity, v, and ion number density, 
n, onboard moments from the electrostatic analyzer (ESA; 
iMcFadden et alJl2008l) . A despun spacecraft coordinate sys- 
tem (DSL) was used, in which z is the spacecraft spin axis. 
The DSL system for ARTEMIS is approximately the same as 
the geocentric solar ecliptic (GSE) system with the sign of the 
>'- and z- axes reversed. The wire boom electric field antennas 
are in the x-y plane and extend a few Debye lengths from the 
spacecraft. Data with the currently most recent calibrations 
(vOl) were used for all instruments. 

For E sc it was found that some extra calibration was needed. 
A least-squares fit, varying the E ytSC offset Oe t and the E sc 
scaling factor F, was performed to minimize the difference 
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Figure 1. Sample power spectra from P2 2010 days 217-230. The dotted 
lines show the range of scales that the spectral index is fitted to. Gradients of 
—5/3 and —3/2 are marked for reference. 

between E y iSC and the y-component of -v x B for each inter- 
val (this technique assumes ideal MHD). Each 6 hr interval 
was corrected using these empirically determined values. The 
mean value of Oe x was found to be -0.17 mV m" 1 for PI 
and -0.23 mV m" 1 for P2; the mean value of F was found to 
be 1.02 for PI and 0.99 for P2. An alternative fit using a B z 
offset instead of an Ey )SC offset was also tried, resulting in B, 
offsets ~ 0.6 nT. The results of this Letter, however, are not 
significantly affected by either of these additional calibration 
methods. 

The velocity and density data was cleaned up by removing 
unphysical spikes and other spurious data that was present. 
This was done by linearly interpolating over data points more 
than 4 standard deviations from the mean in each 6 hr inter- 
val (this process was repeated three times for each interval). 
Any data gaps in the 3 s resolution data were also linearly 
interpolated over to produce time series with consistent 3 s 
resolution. After this process, occasional small spikes, some- 
times seen in all three instruments, remained in the time series 
of some intervals. These are likely due to noise but did not 
affect this analysis (excluding these intervals did not signifi- 
cantly change the results of this Letter). 

The electric field was measured in the frame of the space- 
craft, E sc , and converted into the frame of the mean solar wind 
velocity, E sw , using the Lorentz transformation, 

E sw = E sc + v sw xB, (1) 

where v sw is the mean solar wind velocity relative to the 
spacecraft over each interval (v = v sw + <5v). Since B is a fluc- 
tuating quantity, it was linearly interpolated onto the times of 
E sc so that the transformation could be done for each electric 
field measurement. 

The power spectrum of each component of B and v, of the 
x- and y- components of E sc and E sw , of the magnetic field 
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Figure 2. Histograms of spectral index for the trace of the magnetic field and 
the y-component of the electric field in the spacecraft and mean solar wind 
frames. Values of —5/3 and —3/2 are marked for reference. 

magnitude |B|, and of n was calculated for each of the in- 
tervals. The multitaper m ethod with time-bandwid th product 
AW = 4 and 7 eigentapers dPercival & Waldenl 19931) was used 
(using a standard Fourier transform does not affect the results 
to within errors). Typical power spectra for a longer interval 
(14 days) are shown in Figure Q] the trace of the magnetic 
field spectrum and the y-component of the electric field spec- 
trum in both frames. Since th e solar wind fluct uations are 
anisotropic with k± > k\\ (e.g.. lChen et alJl2010bl) . these are 
measurements of the perpendicular spectrum E(k±) (to mea- 
sure £(^ 11 ) a local field tracking technique would be needed 
(e.g. jHorburv et al.ll2008t IChen et al.ll20lTl) ). 

Each spectral index was determined from the gradient of 
the best-fit line to the power spectrum in log-log space over 
the spacecraft frequency range 1 x 10~ 3 Hz to 2 x 10~ 2 Hz 
(marked as dotted lin es in Figure [T}. Applying Taylor's hy- 
pothesis (iTayloril 1 9 3 8h since the solar wind is super- Alfvenic, 
this range corresponds approximately to scales 18,000 km to 
350,000 km and a perpendicular wavevector range 0.0018 < 
k±pi < 0.036, where p, rj 100 km is the typical ion gyrora- 
dius. This is in the middle of the inertial range and was cho- 
sen because good power laws exist here in all intervals. The 
results of the analysis are described in the next section. 

3. RESULTS 

Histograms of the spectral indices for the magnetic and 
electric fields are shown in Figure [2] The magnetic field trace 
spectral index histogram can be seen to p eak close to -5/3 , 
in agreement with previous results (e.g.. iSmith et al.l 120061 
iTessein et alj 1200 9). The histogram of the y-component of 
the electric field in the spacecraft frame also peaks near -5/3 
but the histogram of the same component in the mean solar 
wind frame peaks closer to -3/2. 
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Table 1 

Mean Spectral Indices 



Field 


Spectral Index 


^ trace 


-1.67 ±0.01 


Vtrace 


-1.50±0.02 




-1.66±0.01 


Ex,sc 


-1.45 ±0.01 


Ey, sw 


-1.40±0.01 


E_x,sw 


-1.39 ±0.01 


|B| 


-1.64±0.01 


n 


-1.63 ±0.02 



The mean spectral indices for each field are given in Ta- 
ble[TJ along with the standard error of the mean cr/y/N, where 
a is the sample standard deviation and TV is the number of in- 
tervals. The mean velocity and density spectral indices were 
calculated from only 117 and 120 of the intervals, respec- 
tively. These are the intervals for which 3 s onboard moment 
data was available and no more than 5% of the data was miss- 
ing. This explains the larger error for these fields. The same 
analysis was also tried with 24 hr intervals (not shown here), 
resulting in a smaller spread of spectral indices but the same 
mean values to within 2 standard errors. 

The fact that the scaling of Ey :SC matches B t race can be shown 
to be due to the Lorentz transformation. In ideal MHD, the 
three fields are related to each other in the mean solar wind 
frame by E sw = — 5v x B and putting this into Equation ([TJ 
gives 

E sc = -(5vxB-v sw xB. (2) 

Since the mean solar wind speed is much larger than the fluc- 
tuations, |v sw | > \5\\, and is mostly in the x (radial) direction, 
Ey iSC is dominated by the magnetic field fluctuations convected 
by the mean solar wind flow and therefore follows their scal- 
ing. The amplitudes of the spectra in Figure [TJ agree with 
this interpretation: the E y>sc spectrum is an order of magni- 
tude larger than the E y syi spectrum, showing that for the y- 
component, the second term on the right-hand side of Equa- 
tion (O is larger than the first. The x-component of Eq. [2] 
does not depend on the radial solar wind velocity, so scaling 
of E XK does not depend only on the scaling of B and indeed 
is different to that of B trace- 

E sw has a scaling closer to that of Vtrace> which can also 
be shown to be consistent with Alfvenic fluctuations in ideal 
MHD. Splitting the magnetic field into a constant mean value 
plus fluctuations, B = Bo + SB, the electric field in the mean 
solar wind frame is given by 

E sw = -(5vx(B + (5B). (3) 

The mean value of |<SB|/|Brj| is between 0.1 and 0.4 for the 
range of scales to which spectral indices were fitted. Since 
at these small scales in the solar wind Bo > SB, the electric 
field spectrum in the mean solar wind frame is dominated 
by the velocity fluctuations, and therefore has a similar scal- 
ing. Similar arguments can be made for the Alfvenic fluctu- 
ations in gyrokinetic theory (Sche kochihin et aT]|2009l) . The 
fact that we observe a spectral index close to -3 /2 in E XiSW 
and £ ( sw also suggests that the perpendicular velocity compo- 
nent has this scaling, whi ch is in agreement with the results of 
Chap man & Hnatl ([2007). The electric field scaling is also in 
agreement with prev ious measurements of the velocity trace 
spect ral index (e.g., Tes sein et al.ll2009t iPodesta & Borovskvl 
l2OT0l) . 



Table 2 

Correlations Between Spectral Indices 



Field 1 


Field 2 


Correlation Coefficient 


^ trace 


E*y,sc 


0.82 


B trace 


Ey } sw 


0.15 


B trace 


Ex.sc 


0.24 


^ trace 


E^sw 


0.08 


Ey.sc 


Ey } SW 


0.17 


Ey } $C 


E x ,sc 


0.24 


Ey^sC 


Ex.sw 


0.14 


Ey.sw 


E X} sc 


0.38 


£y,sw 


Ex,SW 


0.36 


Ex,SC 


E x ,sw 


0.86 



The scaling of the compressive fluctuations (|B| and n) is 
close to -5/3, matching the trace magnetic field spectrum, 
rather than the velocity spectrum. Previous observations (e.g. , 
iMarsch & Tull990tlBellamv et al.l200lllssautier et alj2010h 
could not distinguish between -5/3 and -3/2 in the compres- 
sive fluctuations so this scaling is consistent with those obser- 
vations. Th e compressive fluc tuations are mainly due to the 
slow mode dHowes et al.l 1201 fl) and are sometimes thought 
to be passive to the Alfvenic turbulence. Since their scal- 
ing matches the magnetic field, rather than the velocity, the 
nonlinearity cannot be due solely to passive convection and 
may include nonlinearities with the magnetic field. This sup- 
ports the theori es of compressible reduce d MHD and kinetic 
reduced MHD ( Schekoch ihin et alll2009l) . in which the com- 
pressive fluctuations interact nonlinearly with both the mag- 
netic field and velocity. 

To test the significance of the difference between the mean 
spectral index values in Table [TJ the f-test was applied. This 
is appropriate since the spectral indices appear to be normally 
distributed and are independent measurements. The t value 
for differentiating between the spectral indices of B t race and 
£\>.sc is f = 0.41. This is smaller than the 95% value of 1.96 
for infinite degrees of freedom, showing that there is no sta- 
tistically significant difference between the scaling of these 
two fields. For differentiating between the spectral indices of 
E y>sv/ and E y>sc the t value is t = 15, larger than the 95% value, 
showing that these two fields have significantly different spec- 
tral indices. This confirms that the -5/3 and -3/2 difference 
is a statistically robust result. 

To examine the cause of the spread of spectral index val- 
ues, the correlation between the different spectral indices was 
measured. The linear correlation coefficients, calculated from 
various pairs of sets of the 272 spectral index values of each 
field, are shown in Table [2] It can be seen that the spectral 
indices of most pairs of fields are poorly correlated, having 
correlation coefficients lower than 0.4. This suggests that the 
spread of values is mostly due to random, rather than system- 
atic, variation, although the fact that the correlation coeffi- 
cients are all slightly positive suggests perhaps some small un- 
derlying systematic variation. The exceptions are correlations 
between B t race and E y>sv/ and between E X K and E x>svl , which 
have correlation coefficients larger than 0.8. This is due to 
the reasons discussed above: the E y , sw spectrum is essentially 
a measure of the B t race spectrum because the y-component of 
the last term in Equation (O is large and E XjSC and £ t sw are 
similar because the x-component of the last term in Equation 
(f2]) is not large (since v sw is mostly in the x-direction). 

4. SUMMARY AND CONCLUSIONS 
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We have performed the first survey of electric field data in 
the solar wind to measure the perpendicular spectrum of iner- 
tial range fluctuations. It was found that there is a spread of 
spectral index values but that the average spectral index de- 
pends on the frame of measurement. In the spacecraft frame, 
the j-component of the electric field is primarily due to the 
magnetic field fluctuations being convected past the space- 
craft at the average solar wind speed. It, therefore, has the 
same average spectral index (to within errors) as the mag- 
netic field of —1.66 ±0.01. This is consistent with previous 
single interval electric field measurements in the spacecraft 
frame (|Bale et alJ 12005b ISahraoui et ail 120091) . In the mean 
solar wind frame, the electric field is primarily due to veloc- 
ity fluctuations in a mean magnetic field and has a spectral 
index of — 1.40±0.01, which is closer to the velocity spectral 
index than the magnetic field spectral index, although not the 
same to within errors. The compressive fluctuations (|B| and 
n) were found to have the same spectral index as the magnetic 
field and not the velocity. 

The difference between the scaling of the electric field in 
the spacecraft frame and the mean solar wind frame provides 
independent confirmation of the difference in scaling between 
the velocity and magnetic field. This difference is not ex- 
pected for Alfvenic fluctuations, since 8v is proportional to 
SB in an Alfven wave, and is not predicted by any of the 
current t heories of Alfvenic tu rbul ence (although see recent 
work b vlBoldvrevetal.l (1201 lb and IWang et"aTfil201 11) 1. Re- 
cently, Ro berts! ( 120101) found that further out into the helio- 
sphere, past 5 AU, the velocity spectral index evolves toward 
-5/3 to match the magnetic field. Although an important re- 
sult, this does not explain the difference at 1 AU. Possible rea- 
sons for the difference include the effects of scale-dependent 
alignment, imbalance and residual energy and these will be 
investigated in a future paper. 
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